Light field distribution of general function photonic crystals 
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In this paper, We have presented a new general function photonic crystals (GFPCs), which 
refractive indexes are line functions of space position in two mediums A and B, and obtain new 
results: (1) when the line function of refractive indexes is up or down, the transmissivity can be 
far larger or smaller than 1. (2) when the refractive indexes function increase or decrease along the 
direction of incident light, the light intensity should be magnified or weaken, which can be made 
optical magnifier or attenuator. (3) The GFPCs can be made optical diode when the light positive 
and negative incident the GFPCs. 
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1. Introduction 

Photonic crystals (PC) are a new kind of materials 
which facilitate the control of the light [1-5]. An im- 
portant feature of the photonic crystals is that there 
are allowed and forbidden ranges of frequencies at which 
light propagates in the direction of index periodicity [6- 
9]. Due to the forbidden frequency range, known as pho- 
tonic band gap (PBG) [10-15], which forbids the radia- 
tion propagation in a specific range of frequencies. The 
existence of PBGs will lead to many interesting phenom- 
ena. In the past ten years has been developed an inten- 
sive effort to study and micro-fabricate PBG materials in 
one, two or three dimensions, e.g., modification of spon- 
taneous emission [16-19] and photon localization [20-23]. 
Reduction or suppression of the density of states within 
the band gap facilitates light localization and trapping in 
a bulk material [24-25], as well as the inhibition of spon- 
taneous emission over a broad frequency range. Thus 
numerous applications of photonic crystals have been 
proposed in improving the performance of optoelectronic 
and microwave devices such as high-efficiency semicon- 
ductor lasers, right emitting diodes, wave guides, optical 
filters, high-Q resonators, antennas, frequency-selective 
surface, optical limiters and amplifiers [26-27]. 

In Ref. [28], we have proposed special function pho- 
tonic crystals, which the medium refractive index is the 
function of space position, but the function value of re- 
fractive index is equal at two endpoints of every medium 
A and B. In Ref. [29], we have proposed a new general 
function photonic crystals, i.e., the medium refractive in- 
dex is a arbitrary function of space position, and find the 
transmissivity of one-dimensional GFPCs can be much 
larger or smaller than 1 for different slope linearity re- 
fractive index function. In this paper, we choose two lin- 
earity refractive index functions for two medium A and 
B, and give the light field distribution in the GFPCs. 
We obtain some new results: (1) when the line function 
of refractive indexes is up, the transmissivity can be far 
larger than 1. (2) when the line function of refractive in- 



dexes is down, the transmissivity can be far smaller than 
1. (3) when the two-endpoint values of refractive index 
are equal for mediums A and B, the transmissivity T is 
in the range of and 1 (it becomes conventional PCs ). 
So, the conventional PCs is the special case of the GF- 
PCs. (4) when the refractive indexes function increase 
along the direction of incident light, the light intensity 
should be magnified, which can be made light magnifier. 
(5) when the refractive indexes function decrease along 
the direction of incident light, the light intensity should 
be weaken, which can be made light attenuator. (5) The 
GFPCs can be made light diode. The general function 
photonic crystals can be applied to design more optical 
instruments. 

2. The light motion equation in general func- 
tion photonic crystals 

For the general function photonic crystals, the medium 
refractive index is a periodic function of the space posi- 
tion, which can be written as n(z), n(x, z) and n(x, y, z) 
corresponding to one-dimensional, two-dimensional and 
three-dimensional function photonic crystals. In the fol- 
lowing, we shall deduce the light motion equations of the 
one-dimensional general function photonic crystals, i.e., 
the refractive index function is n = n(z), meanwhile mo- 
tion path is on xz plane. The incident light wave strikes 
plane interface point A, the curves AB and BC are the 
path of incident and reflected light respectively, and they 
are shown in FIG. 1. 

The light motion equation can be obtained by Fermat 
principle, it is 

6 [ n(z)ds = 0. (1) 
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In the two-dimensional transmission space, the line ele- 
ment ds is 

ds = ^{dx) 2 + {dz) 2 = v 7 ! + z 2 dx, (2) 
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FIG. 1: The motion path of light in the medium of refractive 
index n(z). 



where z = then Eq. (1) becomes 

S J n{z)y/l + (z) 2 dx = 
The Eq. (3) change into 
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At the two end points A and f?, their variation is zero, 
i.e., Sz(A) = Sz(B) = 0. For arbitrary variation Sz, the 
Eq. (4) becomes 



dn(z) 
dz 

-n(z\ 



j2_^(£) i 2 (1 + i 2 ) - 

dz 

"i2-i 2 z(l + i 2 ) 



1 + i 2 

simplify Eq. (5), we have 
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The Eq. (6) is light motion equation in one-dimensional 
function photonic crystals. 

3. The transfer matrix of one-dimensional gen- 
eral function photonic crystals 

In this section, we should calculate the transfer matrix 
of one-dimensional general function photonic crystals. In 
fact, there is the reflection and refraction of light at a 
plane surface of two media with different dielectric prop- 
erties. The dynamic properties of the electric field and 
magnetic field are contained in the boundary conditions: 
normal components of D and B are continuous; tangen- 
tial components of E and H are continuous. We consider 
the electric field perpendicular to the plane of incidence, 
and the coordinate system and symbols as shown in FIG. 
2. 

On the two sides of interface I, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, there are 
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FIG. 2: The light transmission and electric magnetic field 
distribution figure in FIG.l medium. 



Eq — Ei — E t i + E' r2 . . 

H =H I = H n cos e\ - H' r2 cos B\. 1 ' 



On the two sides of interface II, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, and give 



En — Ej — E i2 + E r 2 



(8) 



Hi i =H'j = H t2 cos B{ - H r2 cos d\ , 
the electric field E t \ is 

E n = E tw e^ XA+k ^\ z=0 = E no e^ n ^ sin6 ' XA , (9) 
and the electric field E i2 is 

E a = E tw e* k '* XB+k '^\ z = b 

j^j ^gi^n(6)(sin 0*xb+cos 9*b) (10) 

Where x a and xb are x component coordinates corre- 
sponding to point A and point B. We should give the 
relation between Ei 2 and E t \. By integrating the two 
sides of Eq. (6), we can obtain the coordinate compo- 
nent xb of point B 
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FIG. 3: The structure (BA) N of the general function photonii 
crystals. 



where k Q = cot B\ and k z = From Eq. (12), there i: 
n(z) > n(0) sin^j. and the coordinate ig is 



x B = x A + 



f 

Jo 



dz 



V( 1 + fc o)($}) 2 



1 



(14 



where & is the medium thickness of FIG. 1 and FIG. 2. 
By substituting Eqs. (9) and (14)into (10), and using th< 
equality 



we have 
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Substituting Eqs. (16) and (18) into (7) and (8), and 



using H = . I j^nE, we obtain 



where 
M B = 
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The Eq. (20) is the transfer matrix M in the medium of 
FIG. 1 and FIG. 2. By refraction law, we can obtain 
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where no is air refractive index, and n(0) = n(z)\ z -Q. 
Using Eqs. (15) and (21), we can calculate cosOj. 

4. The transmissivity and light field distribu- 
tion of one-dimensional general function photonic 
crystals 
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FIG. 4: The two periods transmission figure of light in general 
function photonic crystals. 



In section 3, we obtain the M matrix of the half pe- 
riod. We know that the conventional photonic crystals 
is constituted by two different refractive index medium, 
and the refractive indexes are not continuous on the in- 
terface of the two mediums. We could devise the one- 
dimensional general function photonic crystals structure 
as follows: in the first half period, the refractive index 
distributing function of medium B is n&(z). and in the 
second half period, the refractive index distributing func- 
tion of medium A is n a (z), corresponding thicknesses are 
b and a, respectively. Their refractive indexes satisfy con- 
dition ribip) i=- n a (0), their structure are shown in FIG. 
3, and FIG. 4. The Eq. (20) is the half period transfer 
matrix of medium B. Obviously, the half period transfer 
matrix of medium A is 
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FIG. 5: The input and output light in the GFPCs. 



cos 



ill 



n 2 a (a) 



sin 



2 9?. 



(26) 



In one period, the transfer matrix M is 
M = M B ■ M A 

cos 6b 
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The form of the GFPCs transfer matrix M is more com- 
plex than the conventional PCs. The angle 9(, d{ , Of 1 
and Of 1 are shown in Fig. 4. The characteristic equation 
of GFPCs is 
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Where N is the period number. With the transfer 
matrix M (Eq. (28)), we can obtain the transmission 
and reflection coefficient t and t, and the transmissivity 
and reflectivity T and R, they are 
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FIG. 6: The line refractive index functions in a period. The 
FIG .6(a) is the up line function of refractive indexes, and 
FIG .6(b) is the down line function of refractive indexes. 



following, we give the electric field distribution of light in 
the one-dimensional GFPCs, and the propagation figure 
of light is shown in FIG. 5. From Eq. (28), we have 
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where di and d 2 are the thickness of first and second 
medium, respectively, E 2 (z) and H 2 (z) are the electric 
field and magnetic field in the second medium, when 
position z changes we can obtain the electric field and 
magnetic field in other medium. The E t and H t are the 
transmission electric and magnetic field. Eq. (33) can be 
written as 

E(z) = A{z)Et N +i+B{z)H tN+ i 
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5. Numerical result 
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In this section, we report our numerical results of trans- 
missivity. We consider refractive indexes of the linearity 
functions in a period, it is 



j^, Eoi and E^ r are incident 



Where rjn = tjn+i 
mil reflection electric field, and E\ = E i + E nr . In the 



n b (z) = n b (0) 



n b (b) - n b (0) 



z, 0<z<b, (36) 



■5 



1.5x10' 



T 

1.0x10 8 - 
5.0x1 7 - 




0.0 I 1 I, I I I ,1 1 ,1 I , II, II, I I I J 1 , 1 I , I I I , 1 I | M , I I 
012345678 

»/co o 



FIG. 7: The relation between transmissivity and frequency 
corresponding to the up line function of refractive indexes 
(FIG .6(a)). 



n a (z)=n a (0)+ riaW a a{ > z, 0<z<a, (37) 

Eqs. (36) and (37) are the line refractive indexes distri- 
bution functions of two half period mediums B and A. 
When the endpoint values n b (0), n b (b), n a (0) and n a (a) 
are all given, the line refractive index functions n b (z) and 
n a (z) are ascertained. The main parameters are: the half 
period thickness b and a, the starting point refractive in- 
dexes n b (0) and n a (0), and end point refractive indexes 
n b (b) and n a (a), the optical thickness of the two medi- 
ums are equal, i.e., n b (0)b = n a (0)a, the incident angle 
9f = 0, the center frequency ujo — 1.215 x 10 15 _ffz, the 
thickness b — 310nm, a — 154nm and the period number 
N = 16, i.e., the structure of GFPCs is (BA) 16 . 

In FIG. 5(a), we take n b {0) = 1.25, n b {b) = 2.18 for 
the medium B, and n a (0) = 2.52, n a (a) = 4.66 for the 
medium A, which are the up line function of refractive 
indexes. In FIG. 5(b), we take n&(0) = 2.18, n b (b) = 1.25 
for the medium B, and n a (0) — 4.66, n a (a) — 2.52 for the 
medium A, which are the down line function of refractive 
indexes. By the refractive indexes function, we can cal- 
culate the transmissivity. With the FIG. 6(a) and 6(b), 
we obtain the transmissivity distribution in FIG. (7) and 
FIG. (8) respectively. From FIG. (7) and FIG. (8), we 
obtain the results: (1) when the line function of refractive 
indexes is up, the transmissivity can be far larger than 
1 (T maximum is 10 8 ). (2) when the line function of 
refractive indexes is down, the transmissivity can be far 
smaller than 1 (T maximum is 1CP 9 ), which are different 
from the transmissivity of conventional PCs (0 < T < 1). 
These results are inexplicable and doubtful even, but 
they are correct. When we take n b (0) = n b (b) = 1.25 
and n a (0) = n a (a) = 2.52, i.e., two-endpoint values of 
refractive index are equal for mediums A and B (it is 
conventional PCs), and by the Eq. (27), (30) and (32), 
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FIG. 8: The relation between transmissivity and frequency 
corresponding to the down line function of refractive indexes 
(FIG .6(b)). 




FIG. 9: The relation between transmissivity and frequency 
for the GFPCs, when the two-endpoint values of refractive 
index are equal 
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FIG. 10: The light positive incident to the GFPCs. 
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FIG. 11: The light negative incident to the GFPCs. 
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FIG. 12: The light distribution of positive incident (FIG. 10) 
in the GFPCs. Figure (a) (b) and (c) are corresponding to 
the incident light frequency uj = loq , 3cj and 5loq ■ 



we obtain the transmissivity distribution in FIG. (9), its 
transmissivity T is in the range of and 1 (it becomes 
conventional PCs transmissivity distribution). So, the 
conventional PCs is the special case of the GFPCs. 

In the following, we shall study the light field distribu- 
tion of the one-dimensional GFPCs for the positive and 
negative incident of light. The positive incident figure is 
shown in FIG. 10, and the negative incident is shown in 
FIG. 11, which relative to the positive incident FIG. 10. 
The refractive indexes line function of positive incident is 
in FIG. 6(a), and the FIG. 6(b) is the refractive indexes 
line function of negative incident corresponding to the 
positive incident FIG. 6(a). From Eq. (35), we can calcu- 
late the electric field distribution in the GFPCs. Figures 
12 and 13 are the electric field distributions of positive 
incident (FIG. 10). In FIG. 12, the frequency of incident 
light uj is the odd times of the center frequency ujq. In 
FIG. 12 (a), (b) and (c), the incident light uj are corre- 
sponding to wo, 3wo and 5wo. When uj = uj and 3uj , the 
light intensity of transmission are enhanced or magnified 
(magnification M = 10 7 ), while uj = 5ujq, the intensity 
of transmission light is weaken. From FIG. 7, we can 
find when uj = oj and 3w the light is in conductance 
band, and the transmission light has been enhanced by 
the GFPCs. In the conventional PCs, the intensity of 
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FIG. 13: The light distribution of positive incident (FIG. 10) 
in the GFPCs. Figure (a) (b) and (c) are corresponding to 
the incident light frequency uj = 2uj , 4cj and 6ujq. 




FIG. 14: The light distribution of negative incident (FIG. 11) 
in the GFPCs. Figure (a) (b) and (c) are corresponding to 
the incident light frequency uj = uio , 3u>o and 5wo • 



transmission light can not be magnified. When uj = 5uj 
the light is in forbidden band, and the transmission light 
has been weaken by the GFPCs. In FIG. 13, the fre- 
quency of incident light uj is the even times of the center 
frequency ujq. In FIG. 13 (a), (b) and (c), the incident 
light uj are corresponding to 2wo, 4wo and Qujq, the inten- 
sities of transmission light have been magnified. Figures 
14 and 15 are the electric field distributions of negative 
incident (FIG. 11). In Figures. 14 and 15, the frequency 
of incident light uj is the odd and even times of the cen- 
ter frequency ujo, and the intensities of transmission light 
have been all weaken. By calculation, we can obtain the 
following results for the GFPCs: (1) For the positive in- 
cident, i.e., the refractive indexes increasing along the 
direction of incident light, the intensity of transmission 
light have been magnified when the incident light fre- 
quency is in conductance band, which can be made light 



7 



0.0 





| (a) £0=6(0, 
V IwvuAA.. _ 


1 1 1 1 1 1 1 1 1 
12 3 4 


? 


(b) to=4co 

lW- i , , , , , ,— 





12 3 4 


■I 


L 





12 3 4 



FIG. 15: The light distribution ol negative incident (FIG. 11) 
in the GFPCs. Figure (a) (b) and (c) are corresponding to 
the incident light frequency uj = 2cj , 4ojo and 6wo- 

magnifier (magnifying multiple 10 7 ).(2) For the negative 
incident, i.e., the refractive indexes decreasing along the 
direction of incident light, the intensity of transmission 
light have been weaken, which can be made light atten- 
uator (magnifying multiple less than 1). (3) The GFPCs 
structure (positive incident FIG. 10 and negative incident 
FIG. 11) is the optical diode, since the light intensity is 
magnified (the light get across) when the light is positive 



incident, and the light intensity has been decreased (the 
light cutoff) when the light is negative incident, which 
transmits light from an input to an output, but not in 
reverse direction. 

6. Conclusion 

In summary, We have theoretically investigated a new 
general function photonic crystals (GFPCs), which re- 
fractive index is an arbitrary function of space position. 
We choose the line refractive index function for two 
mediums A and B. By the calculation, We obtain the 
following results: (1) when the line function of refractive 
indexes is up, the transmissivity can be far larger than 
1. (2) when the line function of refractive indexes is 
down, the transmissivity can be far smaller than 1. (3) 
when the two-endpoint values of refractive index are 
equal for mediums A and B, the transmissivity T is in 
the range of and 1 (it becomes conventional PCs ). So, 
the conventional PCs is the special case of the GFPCs. 
(4) when the refractive indexes function increase along 
the direction of incident light, the light intensity should 
be magnified, which can be made light magnifier. (5) 
when the refractive indexes function decrease along the 
direction of incident light, the light intensity should be 
weaken, which can be made light attenuator. (5) The 
GFPCs can be made light diode. The general function 
photonic crystals can be applied to design more optical 
instruments. 
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